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We study neutrino emission from a normal, ungapped, quark phase in the core of a compact 
star. Neutrino emission from non-interacting quark matter leads to an emissivity that scales as 
e ~ T 7 . We show that the emissivity is enhanced by a combination of Fermi liquid and non- Fermi 
liquid effects. Fermi liquid effects lead to an emissivity that scales as e ~ a s T 6 , as originally shown 
by Iwamoto. We demonstrate that non-Fermi liquid effects further enhance the rate, leading to 
e ~ ct^T 6 log(ra/T) 2 , where m is the electric screening scale and m 2> T under the conditions found 
in compact stars. We show, however, that combined with non-Fermi liquid effects in the specific 
heat the enhancement in the emissivity only leads to a modest reduction in the temperature of the 
star at late times. Our results confirm existing bounds on the presence of ungapped quark matter 
in compact stars. We also discuss neutrino emission from superconducting phases with ungapped 
fermionic excitations. 



I. INTRODUCTION 

Compact stars provide a unique opportunity to study 
cold and dense strongly interacting matter and its rich 
phase structure. An important source of information 
about the structure of compact stars is their cooling be- 
havior 0, E H H @. For the first - 10 5 years after 
the star is born neutrino emission from the bulk is the 
most efficient energy loss mechanism. Since the matter 
in the interior of the star is almost degenerate the ex- 
act nature of the cooling mechanism is very sensitive to 
the structure of the low energy excitations. As a result 
the cooling behavior places important constraints on the 
phase diagram of dense matter. 

Cooling mechanisms are generally grouped into fast 
processes with an emissivity e ~ T 6 and slow processes 
with e ~ T 8 . Slow processes include the modified Urea 
process @ and neutrino bremsstrahlung. Examples for 
fast mechanisms are the direct Urea process 0, neutrino 
emission from pion or kaon condensates 0, U, and neu- 
trino emission from an ungapped quark phase [T(i| . In a 
fully gapped phase of nuclear or quark matter the emis- 
sivity is exponentially small. Near the critical tempera- 
ture for superconductivity the emissivity is dominated by 
pair breaking and recombination with a temperature be- 
havior that is intermediate between fast and slow mech- 
anisms [H03- 

In this paper we wish to focus on the emissiv- 
ity of a possible quark matter phase. At asymptoti- 
cally high baryon density quark matter is in the color- 
superconducting CFL phase ^3|. ^ n this phase all 
quark excitations are gapped and the emissivity is domi- 
nated by exponentially small processes involving massive 
pseudo-Goldstone modes 0, 0, ^| . At densities that 
are relevant to compact stars distortions of the CFL due 
to the non-zero strange quark mass cannot be neglected 
and the phase structure is much more complicated. Most 
of the phases that have been proposed in the literature 
involve quarks with vanishing or very small gaps. Exam- 
ples are the 2SC phase 0, |3 E3 m which one of the 
colors remains ungapped, single flavor spin-one supercon- 



ductors that have very small gaps [2(1 E| , and gapless 
CFL or 2SC phases [HlH. 

In the following we will concentrate on the neutrino 
emissivity of ungapped quark matter. This problem is 
important both in connection with the (almost) gapless 
phases mentioned above and as a benchmark in order 
to exclude the presence of ungapped quark matter in 
compact stars. The emissivity of a normal quark phase 
was first studied by Iwamoto 01 > see a l so 0> HI] • 
Iwamoto showed that the direct Urea process is strongly 
suppressed in non-interacting quark matter, but that 
Fermi liquid corrections in interacting matter lead to a 
fast rate e ~ ot s T & , where a s is the strong coupling con- 
stant. 

It is also known, however, that unscreened color mag- 
netic interactions lead to a breakdown of the Fermi liquid 
description at temperatures T n fi ~ m exp(— 97r/(4a s )), 
where m 2 = Npa s n 2 /tt is the electric screening scale 
p?l |28| . This scale is very small compared to ordinary 
QCD scales, but it is large compared to the temperature 
of neutron stars after the first minute or so. For fi = 500 
MeV we have T n ji ~ 500 keV. It was recently shown that 
non-Fermi liquid effects lead to large corrections of the 
specific heat of degenerate quark matter for T <C T n fi 
|29j . In the present work we compute non- Fermi liquid 
corrections to the emissivity and study the cooling be- 
havior of degenerate ungapped quark matter. 



II. QUARK DISPERSION RELATION 

In a Fermi gas of free quarks the direct Urea process 
is strongly suppressed. As a consequence, the rate is 
very sensitive to modifications of the quark dispersion 
relation. Iwamoto noticed that Fermi liquid corrections 
to the relation between the Fermi energy and the Fermi 
momentum lead to a significant enhancement of the neu- 
trino emissivity ^(j- Non Fermi liquid effects due to 
unscreened transverse gauge boson interactions dramat- 
ically alter the dispersion relation in the vicinity of the 
Fermi surface jH 0, [H EJ HI IH H| and thereby 
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FIG. 1: Feynman diagrams that appear in the matching pro- 
cedure for the forward scattering amplitude to leading order in 
the coupling constant. The upper panel shows the direct and 
exchange terms in QCD. The lower panel shows the collinear 
and contact terms in the effective theory. 



could have a similar effect. In this section we wish to 
study the interplay of Fermi liquid and non-Fermi liq- 
uid effects in dense quark matter and derive the quark 
dispersion relation. 

We are interested in the propagation of quarks in 
the vicinity of the Fermi surface. Since the momentum 
p ~ Vjj, is large, typical soft scatterings cannot change 
the momentum by very much and the velocity is approx- 
imately conserved. An effective field theory for particles 
and holes moving with velocity v = (X,vf) is given by 
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C = ^ v (ivD + 5ti)Tjj v --G^G^+£ HDL 



(i) 



Here, is the covariant derivative and Sfi is a coun- 
terterm for the chemical potential. The effective theory 
contains two types of four fermion operators, correspond- 
ing to forward scattering v + v' — ► v + v' and the BCS 
process v+(~ v) — > v'+(—v') 33]. Since we are interested 
in the normal phase only the forward scattering opera- 
tors are included in equ. (I). The matrix T determines 
the spin, color and flavor structure of the operator and 
w iy-v') is a set of orthogonal polynomials. Equ. QJ also 
contains the hard dense loop effective action 
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which accounts for damping and screening caused by 
particle-hole pairs on the entire Fermi surface. 

The parameters in the effective lagrangian are deter- 
mined by matching QCD Green functions near the Fermi 
surface. In order to determine the dispersion relation we 
have to match forward scattering amplitudes, see Fig.^ 
To leading order in the coupling constant the forward 
scattering amplitude is the sum of a direct and an ex- 
change term. In the effective field theory the direct term 



FIG. 2: Leading order contributions to the fermion self en- 
ergy in the effective theory. The triangle denotes a gluon self 
energy insertion and the square is the four fermion operator 
defined in Fig. 



is reproduced by the collinear interaction while the ex- 
change term has to be matched against a contact term 
|33l |36j . The spin-color-flavor symmetric part is given by 
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with C F = (N 2 - 1)/(2N C ) and all other /? = 0. The 
parameters vf and Sfj, are determined by computing 
the contribution to the fermion dispersion relation from 
states far away from the Fermi surface. This can be done 
most easily using the hard dense loop approximation. We 
find 



Po = \P\ 



which gives 
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These relations can also be derived using the Landau the- 
ory of Fermi liquids. Landau showed that Galilei invari- 
ance implies a relation between the effective interaction 
on the Fermi surface, encoded in the Landau parameters 
/*, and the parameters vf and Sfi. These arguments 
were generalized to the relativistic case by Baym and 
Chin 37]. They show that 



VF_ 



Nf! 



p F = M 1 - 



(6) 



where N = N c Nf/j, 2 /Tr 2 is the density of states on the 
Fermi surface. Using equ. 10 and JBJ gives equ. JSJ. 

Finally, we can study loop corrections in the effective 
theory, see Fig. [21 The collinear loop gives 
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where u> = po — /j,. We recently showed that this re- 
sult does not receive large logarithmic corrections of the 
form g 2n log n (uj) as oj — ► |3J]. The tadpole diagram is 
linearly divergent in the effective field theory. A naive 
estimate can be obtained by cutting the divergence off 
at A ~ Li, This gives E ~ fgNfj, which agrees with the 
result for Sfi given above. We can now summarize the 
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FIG. 3: Kinematics for the quark direct Urea process in the 
limit T ~ E v <C In a free quark gas (upper panel) energy 
momentum conservation forces the quark and electron mo- 
menta to be collinear. If Fermi liquid corrections are taken 
into account (lower panel) the outgoing quark has a non-zero 
transverse momentum l\ ~ a s fig. The dashed dispersion re- 
lations give the HDL result whereas the solid lines show the 
change when non- Fermi liquid corrections are included. These 
lead to a flattening of the dispersion relation in the vicinity 
of the Fermi surface. 



results obtained in this section. The dispersion relation 
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where / is the longitudinal momentum, is given in 

equ. (J7J) and v F ,Sfj, given in equ. (0. 



III. NEUTRINO EMISSIVITY 

The dominant contribution to the emission of neutri- 
nos is given by the quark analogs of /3-dccay (f3) and 
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It is straightforward to introduce weak interactions into 
the effective theory. The charged current interaction is 
given by 
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(ii) 



where cos# c is the Cabbibo angle and gi is related to 
the Fermi coupling by G F /V2 = gl/(8M^). The de- 
pendence on the Cabbibo angle suppresses the processes 
involving the strange quark [lOj. Therefore we will ne- 
glect the neutrino emission of strange quarks. We have 
seen in the previous section that almost collinear gluon 
exchanges can generate large logarithmic corrections to 
the fermion self energy. This raises the question whether 
gluon corrections to the weak interaction vertex have to 
be taken into account. The problem of vertex correc- 
tions was studied in [3ll l3^ | . It was shown that the ver- 
tex receives large logarithmic corrections in the time-like 
regime, but not in the space-like regime. The coefficient 
of the logarithm in the time-like regime is exactly equal 
to the logarithmic term in the fermion self energy. In 
the fermion self energy problem the vertex is space-like 
and does not have to be renormalizcd. However, since 
the wave function renormalization is logarithmically di- 
vergent the effective coupling constant goes to zero near 
the Fermi surface. In the neutrino emission problem the 
vertex is time-like and receives logarithmic corrections. 
However, since the vertex correction is equal to the field 
renormalization the effective coupling constant does not 
change. 

The neutrino emissivity is given by the total energy 
loss due to neutrino emission averaged over the initial 
quark helicities and summed over the final state phase 
space and helicities 
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The weak matrix element for the P and ec processes is given by 

I \ M 0/ed 2 = MG 2 F cos 2 e cP 2 F {vp e ){vp u ), (13) 
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where p e ,Pv are the momenta of the electron and the 
neutrino. Weak processes establish [3 equilibrium in the 
star. In three flavor quark matter with a massive strange 
quark the resulting electron chemical potential is small. 
In the following we shall assume that (T ~ E v ) <C (/i e ~ 
E e ) -C pf- This assumption is appropriate in all cases 
except during the first few seconds of the proto-neutron 
star evolution. 

In this case we can neglect the neutrino momentum 
when applying the energy-momentum conservation rela- 
tion to the matrix element. As a consequence we find 
(v ■ p u ) ~ E v after averaging over the direction of the 
outgoing neutrino. The matrix element is mainly deter- 
mined by the factor (v ■ p e ). To leading order in the 



effective theory the weak decay is exactly collinear and 
(v ■ p e ) = [E e — vple) — up to terms of order 0(T/fj, e ), 
see Fig. If corrections to the dispersion relation are 
taken into account we get 

j2 C (X 

{v-p e ) ~ 0--V F )le + -^T - Sfld-SfJ-u ^ -A*ei ( 14 ) 

Zt e TT 

where l± is the transverse momentum of the electron. 
This result agrees with Iwamoto's result. Non-Fermi liq- 
uid corrections only appear in the phase space integral. 
To leading order in the T//i the sum of the rates for 
electron capture and /3 decay is given by 
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The expression in the square brackets is determined by 
the quark dispersion relation given in equ. ©• Terms 
of 0(a s (a s log(T))™) with n = 0, 1, 2 are independent of 
log(xi) and involve the integral 
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At leading order in T j \x the neutrino emissivity from the 
quark direct Urea process is given by 



457 



G 



Cf a s 



log 



630 " -r-qr-* y ^ 

_(17) 

The first term is the standard result by Iwamoto [lOj , and 
the logarithmic terms are non-Fermi liquid corrections. 
We note that these terms have to be included because at 
very low temperature a s log(T) becomes large compared 
to one. We also note that if the scale inside the loga- 
rithm is on the order of the screening scale, A ~ gfi, then 
a s (/i) log(A/T) stays finite in the limit fj, — > oo at fixed 
T. 



IV. COMPACT STAR COOLING 

In this section we wish to study the impact of non- 
Fermi liquid effects on the cooling history of an isolated 
quark phase. Our aim is not to provide a thorough anal- 
ysis of the cooling behavior of an actual quark or hybrid 
star, but to give a numerical estimate of the size of the 



non-Fermi liquid corrections. The thermal evolution of 
the star is governed by the neutrino emissivity, the spe- 
cific heat and the thermal conductivity. Non- Fermi liquid 
corrections to the specific heat were initially considered 
by Holstein et al. [27| in the case of QED. The calcula- 
tion was recently refined and extended to QCD by Ipp et 
al. They find 
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where the first term is the free gas result and the second 
term is the non-Fermi liquid correction. Ipp et al. also 
determined the scale inside the logarithm as well as frac- 
tional powers of T. From the complete 0(a s ) result we 
find A ~ 0.28 m where m 2 = Nfa s p? /tt is the dynamical 
mass. 

The thermal conductivity of degenerate quark matter 
was studied by Pethick and Heiselberg [33. Their re- 
sult suggests that equilibration is fast and that the quark 
phase is isothermal. In this case the cooling behavior is 
governed by 
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where u is the internal energy, t is time and we have 
assumed that there is no surface emission. Without non- 
Fermi liquid effects we have e ~ T 6 and c v ~ T. In this 
case the temperature scales as T oc l/t~. With loga- 
rithmic corrections included there is no simple analytic 
solution and we have studied equ. (|19fl numerically. 

We take the quark chemical potential to be fi„ = 500 
MeV corresponding to densities ps ~ 6po where po is 
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FIG. 4: Neutrino emissivity e and specific heat c v of quark matter. The dashed line shows the Fermi liquid results and the 
dotted lines show the anomalous corrections. The solid line gives the sum of the two contributions and the gray band shows 
an estimate of the uncertainties. 



nuclear matter saturation density. We note that both c v 
and e are proportional to /x 2 and the main dependence 
of the cooling behavior on fi cancels. We evaluate the 
strong coupling constant using the one loop renormal- 
ization group solution at a scale [i. We take the scale 
parameter to be Aqcd — 250 MeV which gives a s ~ 1 
at fi = 500 MeV. It is clear that the naive use of per- 
turbation theory is in doubt if the coupling is this large. 
In practice we estimate the uncertainty by varying a s 
between 1 and 0.4 which is the value used by Iwamoto 
[To| . We take the scale of non- Fermi liquid effects to be 
A = 0.28m as explained above and assess the uncertainty 
by varying A within a factor of 2. Finally, we took the 
initial temperature to be To = 15 MeV. 

The electron chemical potential is determined by the 
requirements of charge neutrality and /3-equilibrium. In a 
non-interacting quark gas we find /i e ~ m 2 /(App). With 
a strange quark mass m s — 150 MeV this relation gives 
fj, e 11 MeV. This result, however, is very sensitive to 
interactions. To first order in a x the chemical potential 
for a massive strange quark is [3!J 

where E Fs = ^p 2 Fs +m 2 . The important point is that 
the 0(a s m 2 ) term is negative and enhanced by a large 
logarithm \og(p F /m s ). The sign is related to the fact 
that the correlation energy changes sign in going from 
the relativistic to the non-relativistic limit. 

Equation lj20|l implies that the strange quark chemical 
potential can become equal to or even smaller than the 
up quark chemical potential. To leading order in ra 2 /p 2 F 



the electron chemical potential is given by 

Ap F \ 7T \m s ) J 

For the values of the parameters given above this equa- 
tion gives a negative electron chemical potential \i e « — 16 
MeV. In this case the quark phase contains a Fermi sea 
of positrons and the quark direct Urea process is 

u^d + e + + v , d + e + ^u + D. (22) 

The neutrino emissivity is again governed by equ. 1|17|) 
where ^ e has to be replaced by — fi e . We observe that 
despite the large correction to /x e the emissivity is not 
strongly affected. The large variation in fi e when pertur- 
bative corrections are included implies, however, that the 
electron chemical potential is very uncertain. In partic- 
ular, there is a possibility that [i e is much smaller than 
m 2 /(4pi?). If a s fi e < T then the neutrino emissivity is 
no longer proportional to a s [x e T 6 but to T r [25| . 

In the following we shall use the value fi e — 16 MeV cor- 
responding to interacting quarks. In Fig. 0] we show the 
temperature dependence of both the neutrino emissivity 
and the specific heat (solid) compared to the Fermi liquid 
result (dashed) . The gray band shows an estimate of the 
uncertainties which are dominated by the uncertainty in 
the value of the strong coupling. For both c v and e the 
anomalous logarithmic terms (dotted) dominate in the 
relevant temperature range and exceed the Fermi liquid 
result considerably. The cooling behavior is controlled 
by the ratio e/c v . Since e ~ log 2 (T) and c„ ~ log(T) this 
ratio is logarithmically enhanced. However, because the 
temperature at late times scales roughly as the fourth 
root of the numerical coefficient in e/c v this logarithmic 
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FIG. 5: Cooling behavior of ungapped quark matter. We show the temperature Tg in units of 10 9 K as a function of the age 
of the star in years. The dashed line shows the Fermi liquid result whereas the solid line gives the result including non-Fermi 
liquid effects with the estimated uncertainty range. Although the non-Fermi liquid corrections to both the specific heat and 
the neutrino emissivity are significant, there is only a modest reduction in the temperature at late times. However, for both 
cases the cooling of quark matter is considerably faster than the cooling of neutron matter via the modified Urea process given 
by the dotted line. 



enhancement only translates into a modest reduction of 
the temperature. This can be seen in more detail in 
Fig. [5] We observe that compared to the Fermi liquid 
result (dashed) the non-Fermi liquid effects (solid) lead 
to a reduction of the temperature at late times which is 
nearly independent of time. The magnitude of the effect 
is on the order of 20%. For comparison, we also show 
the cooling behavior of normal nuclear matter via the 
modified Urea process n + n^n + p + e~ + i? [fj. We 
have chosen the same density and initial temperature and 
the effective baryon masses given in [f|. We clearly see 
the difference between the fast ~ T 6 quark direct Urea 
process and the slow ~ T 8 modified Urea process. 



V. SUMMARY AND DISCUSSION 

In this work, we have discussed the influence of non- 
Fermi liquid effects on the cooling behavior of compact 
stars containing quark matter in the normal phase. Non- 
Fermi liquid effects lead to a logarithmic enhancement in 
both the neutrino emissivity and the specific heat. The 
net result of these two effects is a mild logarithmic en- 
hancement in the cooling rate. As our rate is even larger 
than the Iwamoto rate we confirm and sharpen earlier 
bounds on the existence of ungapped quark matter in 
neutron stars 0, 0> 0, 0] . More quantitative state- 
ments will require detailed studies of realistic models in 
which the quark core is in contact with a hadronic phase 
or an atmosphere. This is beyond the scope of our inves- 
tigation. 
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e~a 3 T 6 log 2 (T) 


C ~a s ^ 2 Tlog(T) 


CFL 


e ~ r 3/2 exp(-mx/r) 


c v ~T 3 


2SC 


e ~ a a T 6 


c v ~ t?T 


2SC+1SC 


e ~ a s T 6 


c v ~ /i 2 T 


[ISC] 3 


e ~ a 3 T 6 


c v ~ t?T 


gCFL 


6 ~ q s T 5 ' 5 


c v ~ /iW 


gCFLK 


e ~ a a T 8 


c v ~ t?T 


g2SC 


e ~ a 3 T 6 


Cv ~ fi 2 T 



TABLE I: Dominant cooling rates in different phases of quark 
matter. This table is a summary of the discussion in Sect.lVI 
In the entry [ISC] 3 we have assumed that ungapped modes 
exist and the Urea process is allowed. This is not true in the 
CSL phase. In the gCFL phase we have used the estimate 
of Kouvaris et al. and in the gCFLK phase we have assumed 
that only slow processes such as neutrino bremsstrahlung are 
allowed. In the g2SC phase we have assumed that the Fermi 
liquid direct Urea rate can be used. In both the gCFL and 
g2SC phase we have ignored possible instabilities signaled by 
a negative current-current correlation function. 



Let us now discuss the importance of non-Fermi liquid 
effects in partially gapped color superconducting phases, 
see Table H] The simplest case is the 2SC phase. The 
2SC phase can arise when the difference between the 
Fermi momenta of the strange quark and the up and 
down quarks is too large for strange-non-strange pairing 
to occur 0, ^(| . There also regions in the phase dia- 
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gram where us (dSC) or ds (uSC) pairing might occur 
H3) HE E3- The 2SC phase is characterized by a par- 
tial Higgs mechanism. Color 5(7(3) is broken according 
to SU(3) — > SU(2) and five out of eight gluons acquire 
a mass. The gapless fermions of the third color inter- 
act via screened gauge bosons. As a consequence, there 
are no non-Fermi liquid corrections in either the specific 
heat or the neutrino emissivity and the standard result 
of Iwamoto applies [l0( ■ 

If the 2SC phase occurs in quark matter with three 
flavors then the unpaired quark flavor can form a spin 
one condensate. Spin one condensates are interesting 
from the point of view of cooling because the typical gaps 
are much smaller than the spin zero gap. In particular, 
we can have T c ~ Tnfi and non-Fermi liquid effects in 
the specific heat are important. Whether or not gapless 
modes occur below T c depends on the exact nature of 
the spin one condensate pol |2]| . In a phase with CSL 
pairing in both LL and LR chirality channels the order 
parameter is isotropic and no gapless modes exist. In all 
other phases not all fermions are gapped and the order 
parameter has nodes on Fermi surface. If spin one pairing 
takes place in all flavor sectors (referred to as [ISC] 3 in 
Table[fl these modes will dominate the cooling behavior. 

The cooling behavior of the recently proposed gapless 
CFL and 2SC j2^| phases is a more difficult question. 
In the gCFL phase there are gapless charged and neutral 
modes and the direct Urea process is possible. Alford 
et al. [5(j have argued that the (almost) quadratic dis- 
persion relation of one of the gapless modes leads to an 
enhancement of the cooling rate by a factor \/T . A sim- 
ilar enhancement occurs in the specific heat |22|. Kry- 
jevski and Schafer studied quark modes in a kaon con- 
densed CFL phase j^Hl- They find at most one gapless 
mode. As a consequence, the dominant cooling mecha- 
nism is expected to be a slow process such as neutrino 
bremsstrahlung. However, there are several light modes 



that will contribute to neutrino transport at MeV tem- 
peratures. In the g2SC phase there are gapless up and 
down quark excitations and the direct Urea process is 
possible. 

We should stress that our discussion of neutrino emis- 
sion from the g2SC and gCFL phase ignores potential 
instabilities related to a negative current-current corre- 
lation function 15 2t 153 , |5J] . Wu and Yip as well as Gian- 
nakis and Ren [52,|55| have argued that these instabilities 
indicate that the correct ground state is an inhomoge- 
neous superconductor of the type discussed by Larkin, 
Ovchninkov, Fulde and Ferrell [il HI ill . The LOFF 
phase also contains gapless fermions and the calculation 
of the neutrino emissivity from this phase in an impor- 
tant problem for future studies. 

Note added in proof: 
Since the first version of this work the scale in the log- 
arithm of equ. has been computed by Gerhold and 
Rebhan and the result As ~ 8m turned out to be 
substantially larger than the corresponding scale in the 
specific heat A c i=s 0.28m we have employed as an esti- 
mate for the scale in the emissivity. This leads quantita- 
tively to a further enhancement of the neutrino emissivity 
but does not change the cooling behavior significantly. 
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